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Abstract: We use a recently proposed non-perturbative model, based on an effective 
strong coupling constant and free from tunable parameters, to study c- flavoured hadron 
production in e~^e~ annihilation. Charm-quark production is described in the framework 
of perturbative fragmentation functions, with NLO coefficient functions, NLL non-singlet 
DGLAP evolution and NNLL large-x resummation. We model hadronization effects by 
means of the effective coupling constant in the NNLO approximation and compare our 
results with experimental data taken at the pole and at the T(4S') resonance. We find 
that, within the experimental and theoretical uncertainties, our model is able to give a 
reasonable description of D*^-meson spectra from ALEPH for x < l — A/nic- More serious 
discrepancies are instead present when comparing with D and D* data from BELLE and 
CLEO in x-space. Within the errors, our model is nonetheless capable of reproducing the 
first ten Mellin moments of all considered data sets. However, the fairly large theoretical 
uncertainties call for a full NNLO/NNLL analysis. 
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1. Introduction 

The hadronization of partons into hadrons, for the time being, cannot be calculated from 
first principles, but it is usually described in terms of phenomenological models, such as 
the Kartvelishvili [|| or Peterson non-perturbative fragmentation functions, containing 
few parameters which need to be tuned to experimental data. It was recently proposed 
m, ^, however, a non-perturbative model, based on the work in refs. [^, |6|, including power 
corrections via an effective strong coupling constant, which does not exhibit the Landau 
pole any longer and includes absorptive effects due to gluon branching. The interesting 
feature of such a model is that it does not contain any extra free parameter to be fitted to 
the data, besides the ones entering in the parton-level calculation. In Q, such a model was 
used in the framework of S-meson decays and it was found good agreement with the data 
on the photon spectrum and on the hadron-mass distribution in radiative and semileptonic 
decays, respectively. In the effective coupling was employed in the framework of bottom- 
quark fragmentation and, within the theoretical uncertainties, a reasonable fit of LEP and 
SLD data on B-hadron spectra was obtained in both x and moment spaces. 

Although the results in [^, |^ are encouraging, it is nonetheless mandatory to con- 
sider more data and observables to validate the effective-coupling model. In this paper, 
we consider charm-quark fragmentation in e'^e~ processes and investigate how our non- 
perturbative model fares against D-meson data from LEP and S-factories. In fact, charm 
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production involves pretty different scales with respect to 6-quark fragmentation, and there- 
fore the comparison with D-hadron spectra should help to shed light on our model. Con- 
sidering charm production at the pole and at the T(4S') resonance, furthermore, is also 
interesting to understand how our model behaves when the process hard scale changes. 
Perturbative charm production will be described in the framework of perturbative frag- 
mentation functions Q , using the same approximations carried out in , and the effective 
coupling constant will be our only source of non-perturbative power corrections. 

The plan of the present paper is the following. In section 2 we shall review the main 
points of the parton-level computation, based on the perturbative fragmentation formalism, 
and including large-x resummation in both coefficient function and initial condition of the 
perturbative fragmentation function. In section 3 we shall discuss the effective coupling 
constant and the inclusion of non-perturbative corrections to charm-quark fragmentation. 
In section 4 we shall compare the results with charmed-meson spectra from LEP and B- 
factories in x-space, whereas we present our analysis in Mellin moment space in section 5. 
We shall finally summarize our main results in section 6. 



2. Charm-quark production 

In this section we shall discuss our calculation for charm-quark production. For the sake 
of consistency, and given the tight relation between perturbative calculation and non- 
perturbative corrections, our computation will be carried out along the lines of ref. Q. 
Therefore, we shall just point out the main issues involved in the calculation and refer to 
Q for further details. 

2.1 Perturbative fragmentation functions 

We consider cc-pair production in e^e~ annihilation at next-to-leading order (NLO) in the 
strong coupling constant as'- 

e+e- ^ P(Q) ^ c(p,)c(pe-) {.9{Vg)) (2.1) 
and define the charm-quark energy fraction: 



2.2 



In the following, we shall consider charm production at LEP, where P \s a. boson and 
Q = mzi as well as c-quark fragmentation at the T(4S') resonance, i.e. Q = ?7i-f(4s) 
and the cc pair coming from the decay of a virtual photon (P = 7*). In principle, in 
e+e~ annihilation charm quarks can also come from other processes, such as the decay of 
bottomed hadrons produced via Z^{T{AS)) — > bb. However, as we shall discuss in detail in 
section 4, our analysis will only deal with direct cc production. 

The perturbative fragmentation approach 0, up to power corrections, factorizes the 
energy distribution of a heavy quark, the charm quark in our case, as the convolution of a 
coefficient function, associated with the emission off a massless parton, and a perturbative 
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fragmentation function, expressing the transition of the hght parton into a heavy quark. 
This way, the c-quark spectrum reads: 



+ 0{{mJQr). (2.3) 



1 d&i 



--7^(^,Q,/^i?,AtF) 
a az 



In eq. p > 1, ddi/dz is the differential cross section for the production of a massless 

parton i after subtracting the cohinear singularity in the MS factorization scheme; ixr 
and lip are the renormalization and factorization scales; a is the NLO e^e~ — > qq{g) 
cross section. Hereafter, we shall neglect charm production via g ^ cc splitting. In 
fact, we can anticipate that, when comparing with data, secondary cc production will 
be either subtracted off the sample which we shall analyse or negligible at the centre-of- 
mass considered energy. This implies that i = c in eq. (^) and is the perturbative 
fragmentation function expressing the fragmentation of a massless c into a massive c. The 
NLO MS coefficient function for e~^e~ — > qq processes can be found in Q. 

The perturbative fragmentation function follows the DGLAP evolution equations ||9|, 
p!o[| ; its value at a any scale fip can be obtained once an initial condition at //qf is given. 
In 1^] the initial condition D™^{x, hqf ,mc) was calculated in the NLO approximation and 
its process-independence was established on more general grounds in |jll|]. It is given by: 

,2 N ,,2 _2n _ , as(MoJ?)Ci^ \l + x^ A„/^OF 



Z)™(x,a5(/ioi?),/xoF>"iJ = 5{l-x) + ' 



\ — X \ ra^ 



\-at3UL-2\n{l-x)-l 



2tt 

(2.4) 

As discussed in [Q], solving the DGLAP equations for an evolution from //qf to fip, 
with a NLO kernel, allows one to resum leading (LL) In^ {iJ.'jp / iJ,Qp) and next-to-leading 
(NLL) a'g ln'^~^ {iJ?p / fj,Qp) logarithms. Setting ^uqf — w-c and //^ ~ Q, one resums the large 
ln{Q'^/mf.) appearing in the massive NLO spectrum 0. The resummation of such mass 
logarithms is usually called collinear resummation. For the sake of working in the same 
perturbative framework as in Q, in the following we shall consider NLO coefficient functions 
and initial condition, along with NLL non-singlet DGLAP evolution. However, one could 
go beyond such a level of accuracy and include NNLO corrections to the coefficient function 
p!^, 13|, initial condition |14] and to the non-singlet splitting functions [15| entering in the 



kernel of the DGLAP equations. The gluon-initiated contribution to the initial condition, 
necessary to possibly extend the analysis to the singlet sector, was calculated in Q and 



\1§ to NLO and NNLO, respectively. 



2.2 Large-x resummation 

Both coefficient function ||7| and initial condition (2.4) contain terms, ~ 1/(1 — x)+ and 



~ [ln(l — x)/(l — x)] + , enhanced when x approaches 1, which corresponds to soft- or 
collinear-gluon radiation. One needs to resum such contributions to all orders to im- 
prove the perturbative prediction (threshold resummation). As in we shall implement 
threshold resummation, which is process-dependent in the coefficient function and process- 
independent in the initial condition |^], in the next-to-next-to- leading logarithmic (NNLL) 
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approximation, following the general method of 18 1. Large-x resummation is typically 
performed in Mellin moment-space, where the Mellin transform of the differential cross 
section reads: 



AT-i 1 da 



(TN = dx x'^~' — - . (2.5) 
a dx 







In iV-space, the enhanced contributions ~ as/(l— x)+ and as'[ln(l— x)/(l— correspond 
to single (~ as In A^) and double (~ as In^ N) logarithms of the Mellin variable A'^. The 
resummed coefficient function is given by the following generalized exponential function 



[as{^%),fi%,lil,Q^] = exp{G^S^ [as{^%), ^^l,Q']} , (2.6) 



where 



G^P[asif.l),f.l,^hQ'] = / dz'-^\ i ^^M^^)] 

JO 

+ B[as{Q\l-z))]]. (2.7) 

The exponent [as'(/i^), /i^, Q^] resums the large logarithms of the Mellin variable; 
in the NNLL approximation, one keeps in the exponent terms ~ a^ln"^^A^ (LL)) ~ 
agln"iV (NLL) and ~ agln""^iV (NNLL). As in the integration variables are z = 
1—Xg, Xg being the gluon energy fraction, and /c^ = [pc+Pg)"^ z) . In soft approximation, 
z ~ x; for small- angle radiation k'^ ~ k\, the gluon transverse momentum with respect to 
the c. 



In (2^), function A[as) resums soft and collinear radiation, while B[as) includes 



all-order collinear and hard emissions. They can be expanded as a series in as as: 

oo 

•*M = E(v)"'""'- P- 



vr 

n=l 



n=l 

In the NLL approximation, one needs to include the first two coefficients of A[as) and the 
first of B^as)] to NNLL accuracy, 

^(3) and 5(2) 

are also needed. The coefficients A^^\ A^"^^ 
and can be found in ||l^; more recent is the calculation of the NNLL contributions 



^(3) H and ||. 

Likewise, the threshold-resummed initial condition reads [p] ]: 

aS?^ [a5(/^o/?)> /^o/?> /^of> "^c] = exp \g^P [a5(/^o/?)> /^of> "^c] } > (2-10) 

where 



Cjv [«5(^D)/^oi?)/^OF>"^c] = / dz— I T2"^["S'(^^)] 

JO I — Z (^7m2(i_2)2 /C 

+ D[asK(l-z)2)]}, (2.11) 



- 4 - 



with /c^ and z defined as in (2/7). To NNLL accuracy, we need J^^^ ^ A^'^^ and A^^^ and the 
first two coefficients of 

oo 
n=l 

namely D^^^ and D^'^^ . Function D{as), called H{as) in |11|, is characteristic of the frag- 
mentation of heavy quarks and resums soft and large-angle radiation. Its 0{as) coefficient 
can be found in [11|, while D*-^^ can be read from the formulas in pTI, 22, |14|. In any case, 



all relevant NNLL threshold-resummation coefficients are reported in |^]. 

In the phenomenological analysis of the inclusion of NNLL effects, and especially 
the contribution ~ a^^A^^^^ to function A{as), turned out to be necessary to reproduce 
the 6-fragmentation data. In fact, as we shall point out in the next section, when using 
the effective coupling constant we need to redefine the threshold-resummation coefficients 
from the third order on. This way, it turns out that A'^'^^ gets enhanced. The inclusion of 
NNLL terms in the resummed exponents shifted the S-hadron spectrum towards lower xb 
values and played a crucial role to obtain a reasonable description of LEP and SLD data 
(see figure 4 in ref. 

As in 1^ §1, the Mellin transforms of our resummed expressions will be performed 
exactly and not according to the step-function approximation, which was instead employed 
in the resummations carried out in refs. [jll], |l^ ^. In fact, as we shall discuss later, we will 
model non-perturbative effects to charm fragmentation by means of an effective coupling 
constant and it was found in ||2^ that the step-function approximation would suppress most 
power corrections included in the physical observables via the analytic coupling. In any 
case, as thoroughly detailed in Q, the issue of the power corrections which are transferred 
to the cross section by the effective coupling, and whether it is a better approximation 
performing the Mellin transforms in an exact or approximated way is currently an open 



issue and we cannot draw any firm conclusion. A careful analysis, along the lines of |27|, 
will be anyway very welcome to clarify this point. For the time being, the exactness of the 
Mellin transforms should be seen as part of our non-perturbative model. We just point 



out that, unlike refs. |11, 18 1, where only logarithms of N are resummed, in our approach 
even some constants and power-suppressed 0{1/N) terms are included in the exponents 
( |2.6| ) and ( p. 10 ) thanks to the exact Mellin transforms. This implies that, any time we 



improve the accuracy of the large-x resummation, e.g., from NLL to NNLL, we include in 
the resummed exponent not only subleading logarithms of A^, but also constants and power 
corrections. This tight relation between perturbative and non-perturbative corrections is 
indeed a peculiar feature of our effective-coupling model. 

As in IQ], the resummed results are matched to the exact NLO coefficient function and 
initial condition. A difference with respect to the approach followed in ref. is that we 
implement the so-called In /^-matching |2^, corresponding to matching the logarithms 
of resummed and NLO expressions. We briefly review this matching strategy and how it 



""^In [|ll| , [L8[ , the longitudinal-momentum integration is done after performing the replacement —^ 
1 — 0^1 — 2: — ° ^ , which is a correct approximation to NLL accuracy. Beyond NLL, it can be generalized 
following the prescription presented in |25[. 
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compares with the standard method implemented in Q. Referring, e.g., to the coefficient 
function, matched to the exact NLO one, it can be written (see eq. (4.2) in ref. Q) as: 

+ d'^P[as{fil),f^%,fil,Q^]. (2.13) 



In ( 2.13| ), is the resummed coefficient function, given in eq. (|2] 

i^(^) [asi^^l),^J.ji,^J.l,Q'] =l + asifil)Qifil,Q') (2.14) 

is a hard factor, including the constant terms which are present in the NLO coefficient 

(c) 

function but are not resummed in A^y ' , 

SP [asipl),fJ.%f^l,Q^] = asifil)Y{f,l,Q^) (2.15) 

is a remainder function, collecting the left-over NLO terms, suppressed at large N. The 
explicit expression for functions Q{n'^,Q'^) and y(^|,,Q^) can be read from the formulas 
in Q. A similar expression holds for the resummed initial condition matched to the NLO 
result (see eq. (5.17) in ref. Q). 

According to the In i?-matching, functions K^'-"^ and d^'-'^ are to be replaced by expo- 
nential functions of their 0{as) terms and eq.( p.l3 ) should read: 



X eMasif^R)yif^F,Q% (2.16) 



From eq. ( 2.16| ), one can easily check that the logarithms of NLO and resummed functions 

v2 - 



are actually matched. In particular, eq. ( |2.16| ) differs from (2.13) only by terms of 0{a'gj 



or higher, but it is smoother at small and large values of (x), thanks to the exponential 
functions in eq. ( 2.16 ). It was in fact pointed out in Q that, since the remainder function 



contains terms ~ Inx and ~ ln(l — x), the physical differential cross sections exhibit 
oscillating behaviour near x ~ and x ~ 1. Exponentiating the 0(05-) contributions to 
the remainder function should therefore improve the prediction for small and large values of 
X. The In /^-matching prescription will be adopted in the following even for the resummed 
initial condition of the perturbative fragmentation function. 

The c-quark spectrum will finally read in A-space as follows: 

c'n [as(/^L)> "s(m|), /^L' iAf^ lA^^ral, Q'^]= C^*^ [as(/i|), /x^^, Q^] 

xEN[as{tilF),as{lA^)\ (2.17) 

^, T-,im,res \ / 2 \ 2 2 21 



In eq. ( p. 17 ), -Eat [a5(/io_F)) ck5(/^f)] is the DGLAP operator for an evolution between the 



scales ^oF and ^p- Throughout this paper, we shall implement E'at [^^(/Uq^), as(/i|.)] in 
the non-singlet approximation; its explicit expression can be found, e.g., in ref. 0. 
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3. Effective coupling constant 



We shall include non-perturbative corrections to charm fragmentation using, as in Q, 
a model, based on an extension of refs. [^, |6|, which includes power corrections via an 
effective strong coupling constant, and does not introduce any further parameter to be 
tuned to experimental data. We review below the main points of our model. 

As discussed in ref. [28|, in resummed calculations the momentum-independent cou- 
pling constant is replaced by the following integral over the discontinuity of the gluon 
propagator: 

«s - 7^ ds Disc, (3.1) 

where k'^ is the gluon transverse momentum relative to the emitter, defined, e.g., as in 
eq. ( p.?! ). In eq. ( |3.lD the discontinuity is given by: 

DisCsF(s) = lim [F(s + ie) - F(s - ie)] . (3.2) 
At LO, e.g., as{—s) reads: 

/3o[ln(|s|/A^) - «7rG(s)] 

where /?o = (33 — 2nj)/(127r) is the first-order term of the QCD /3-function, Uf is number 
of active flavours, and A is the QCD scale, e.g., in the MS renormalization scheme. 

The integral ( |3.1| ) is usually carried out neglecting the imaginary part, ~ in, in the 
denominator of as{—s), i.e. assuming 

ln-p>7r (3.4) 

in eq. ( |3.3| ). The approximation ( |3.4| ) allows one to avoid the Landau pole, so that the 
integral (^]^) turns out to be roughly equal to the strong coupling constant evaluated at 
the upper integration limit: 

— ds Disc, ^^izfl ^ asik^). (3.5) 

In fact, resummed formulas typically use the transverse momentum fc^ as the scale of the 
strong coupling constant p^ . 

As in we shall follow a different approach and avoid the Landau pole by using in 
eq. ( p. ID a regularized coupling constant as, defined as follows [^, ^: 

1 ds 

c^sik') = — ^Disc.«s(-«). (3.6) 



Inserting in (3^) the LO expression (|3.3| ) and performing the integration, we obtain: 

(3.7) 



as,Lo{k ) = -r 



ln(fcVA2) - A2 



-7- 



If we compare eq. (|3.7l) with the LO standard couphng, i.e. 



we learn that in eq. (3/7) a power-suppressed term, relevant at small k'^, has subtracted 
off the Landau pole k'^ = h? , which is instead present in (3^). At large k"^, as{k'^) 



is nonetheless still roughly equal to as{k'^)- Such results can be generalized to higher 
accuracy levels, using the two- and three-loop beta function, as done in |Q. 

The effective coupling constant 05 (/c^) will be still defined as in eq.(^), but using the 
analytic coupling in the integrand function: 



2tt 



asik'') = — ds Disc, ^^L^. (3.9) 







Using the LO result (3.7), we can perform the integral ( |3.9[ ) and obtain our LO effective 
coupling constant: 

"ln(A;VA2) 



/,2^ 1 I 1 1 

as,Lo[k ) = -r "{t; arctan 



vr 



(3.10) 



The NLO and NNLO expressions of as{k'^) can be found in It is straightforward to 
show that eq. ( p. 101 ), as well as its higher-order generalizations, is free from the Landau 
pole and includes power-suppressed contributions at small momenta. Also, as discussed 
in [|^], eq. ( |3.9D accounts for absorptive effects due to gluon branching, since we are not 
neglecting any longer the imaginary part in the denominator of as{—s). 

In principle, both analytic coupling constants (|3.6|) and (|3.9|) are possible candidates 
to model non-perturbative corrections ^. However, as debated in [Q], it is only ( p.S| ) which 
gives an acceptable description of 5-fragmentation data and we shall therefore stick to 
as{k'^) to model power corrections to charm fragmentation as well. 

The relation between effective and standard coupling constant for ln{k'^ /A'^) S> vr reads: 

as{k') = as{k') - M)!a|(fc2) + 0(a|). (3.11) 

From eq. ( |3.11| ) we learn that at high energy the difference between as{k'^) and as{k'^) starts 
from 0{a^). Moreover, eq. ( 3.11| ) dictates that, when employing the effective coupling 
constant, we will have to redefine the soft-resummation coefficients from order a| on. As 
anticipated in subsection |2.2| , the NNLL coefficient ^4^^^ of the 0{a^) term of function 
A{as), entering in Eqs. ( |2.7D and ( 2.11| ), will get enhanced according to: 

^(3) ^ i(3) ^ ^(3) + M)!^(l)^ (3_12) 

The other assumptions contained in our model are also detailed in ref. Q and we do not 
report them here for the sake of brevity. We just point out that, when dealing with higher 



^In the literature one usually refers to Qs(fc^) and Qs(fc^) as efltective space- and time-like coupling 

constants, respectively. 
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orders of asik"^), we shall adopt the so-called 'power-expansion' choice, which implies that 
we shall evaluate the powers a^(/c^) after computing the integral over the discontinuity: 

n 

(3.13) 

On the contrary, the original proposal in |^ consisted in calculating the discontinuity of 
ag{—s) before integrating over s ('non power-expansion' choice). As discussed in Q, the 
non-power expansion prescription would yield a rather poor description of 6-fragmentation 
data. 

The purpose of the present paper is indeed to push the effective-coupling model to 
lower energies and compare its predictions with data on c-flavoured hadron production. 
For a consistent comparison with the results obtained in the framework of i?-hadron pro- 
duction and decay, throughout this paper, as in |3|, we shall use dsik'^) evaluated to 
three-loop accuracy, everywhere in our calculation, i.e. in both coefficient function and 
perturbative fragmentation function. Hereafter, the effective coupling constant (^]^) will 
be our only source of non-perturbative corrections and we shall not introduce any further 
non-perturbative fragmentation function. 

It was pointed out in |Q that power-correction effects in the initial condition of the 
perturbative fragmentation function are more relevant than in the coefficient function. 
The typical c- fragmentation scales at which the coupling constant is evaluated are, in fact, 
C = QVl ~ X in the coefficient function and S = rn-c(l — x) in the initial condition. C and 
S are the integration limits in the resummed exponents as well as the arguments of as in 
functions B [as(C^)] and D [05(5"^)], appearing in the large-x resummation expressions 
( |2.7| ) and ( |2.11| ). If we calculate C and S for Q = mz and x = 0.5, where, as will be 
shown in the next section, the Z)-meson spectrum in e~^e~ annihilation is roughly peaked, 
we shall get C ~ 46 GeV, S ~ 0.9 GeV, dsiC^) ^ 0.13 and ds{S^) ^ 0.35. Therefore, 
non-perturbative corrections are more important in the initial condition, depending on 
S, than in the coefficient function. Comparing now the values of 0:5 at the charm- and 
bottom-mass scales, we find that d5(m^) ~ 0.3 is appreciably higher than Q5(m^) ~ 0.2. 
However, it is interesting to notice that the scales S and C, and hence 0:5(5^) and dsiC^), 
are roughly the same for bottom and charm production if evaluated at the maxima of the 
respective spectra at LEP, i.e. x=0.5 for D- and x = 0.8 for i?-hadron energy distributions. 

Before closing this section, we would like to stress that, in its current formulation, 
our parameter-free model works in the same fashion for B as well D mesons, up to the 
replacement rrib rric- Also, our model does not distinguish among baryons and mesons, 
spin-1 and spin-0, charged and neutral hadrons. It was therefore argued in |^] that possible 
extensions of our model may consist in including a correcting term, so that 

ds{k^) ^ ds{k^) + 5ds{k^) , (3.14) 

where ds{k'^) is still the effective coupling discussed above, and (5qs(A:^) may depend, e.g., 
on whether we have baryons or mesons, -B's or D's, and so on. The analysis which we shall 
undertake herafter should therefore be helpful to establish, for the time being, whether the 
contribution 5ds{k'^) is mandatory or not. 



d^sik^ 



2^ 



ds Disc. 



as{-s) 
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4. Results in x-space 



In this section we compare our results in x-space with experimental data on c-flavoured 
hadron production in e~^e~ annihilation. Hadronization effects will be accounted for by 
employing the analytic coupling constant ( |3.9| ) at NNLO. Whenever we use asik"^) instead 
of the standard as{k'^), the charm-quark energy fraction will be replaced by its hadron-level 
counterpart: 

with p£) being the momentum of a D-hadron. The D spectrum in moment space will be 
written in a form analogous to eq.( p.l7 ), up to as ^ as- 



X Dj^^'"^" [as{^J.lR), f^R, f4F,ml] . (4.2) 
The x-space result is then recovered by performing an inverse Mellin transform: 

/'74-20O J AT 

(D) I 2 2 2 2 2 r^2\ / -A^ {D) , 2 2 2 2 2 r^2\ 

a"- > [xd; f^R,fi0R^fJ-0F,f^F,m^,Q ) = / IT^^D ^iV (/"R' /^Oi?' /^0F> Mf> "^c> Q )> 

(4.3) 

where 7 is a positive constant. As discussed in since the effective as{k'^) does not 
exhibit the Landau pole any longer, we do not need any prescription, such as the well- 



known minimal prescription |29|, to avoid the Landau pole in the integration ( |4.3| ). The 
integral will be performed in a numerical way, along the lines of Q; it was checked that 
the results are stable when varying the integration contour, i.e. the constant 7. 



As in ref. [30|, we shall consider LEP data from the ALEPH collaboration taken 
at the pole, and data from the CLEO and BELLE |33] experiments, at the T(4S') 



resonance. We shall investigate neutral as well as charged D and D* mesons; in fact, we 
just pointed out that our model does not distinguish the hadron electric charge or spin. 

As discussed in [^], electromagnetic initial-state radiation (ISR) effects can modify 
the shape of charmed-meson spectra. Such effects are important especially at i?-factories, 
where the emission of photons from the e~^e~ pair, whose rate is ~ aln{Q'^ /ml), me being 
the electron mass, may significantly decrease the energy in the centre-of-mass system. The 
CLEO and BELLE data did not account for such effects, which were instead implemented 



in the analysis |3C]. In the following, we shall compare with data corrected for ISR effects: 
a discussion on the impact of such contributions on D-spectra in x- and A^-spaces can be 
found in |^^. Such effects were also implemented to correct the ALEPH data, but it was 
understood that at the pole they are quite negligible. 

The non-perturbative model based on the effective coupling constant ( |3.9D does not 
have any free parameter to be tuned to the data which we shall consider. We shall nonethe- 
less vary the parameters entering in the perturbative calculation in such a way to give an 
estimate of the theoretical uncertainty on our prediction. We change each quantity sepa- 
rately, keeping the others to their default values, in such a way to avoid too many runs. 

Following [Q], the default values of our perturbative parameters will be fiR = = Q 
and fiQR = = mc, where fiR and hf are the renormalization and factorization scales 
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in the coefficient function, and hqr and /iQF in the initial condition of the perturbative 
fragmentation function. The hard scale will be Q = mz or mx(4s) at LEP or i?-factories, 
with mz = 91.19 GeV and mx(45) = 10.58 GeV. We shall vary fiR and hf between Q/2 and 
2(5, fJ-oR and between mc/2 and 2mc. As in we shall let as{m'^) run in the range 
0.117 < as(m|) < 0.121, using as{m?z) ~ 0-119 our default value. The corresponding 
variation range of the effective coupling constant is 0.115 < asi'f^'z) ^ 0.119. For the 
purpose of rric, as thoroughly discussed in using the pole or the MS heavy-quark mass 
definition in the initial condition is equivalent for calculations relying on the NLO/NLL 
approximation. However, in the NNLL large-x resummation of the initial condition, and 
in particular in the definition of the coefficient D^'^^ in eq. ( 2.11 ), we are employing results 



of the NNLO computation in |14], which uses the heavy-quark pole mass. Hence, we 
should use the charm pole mass as well. Nonetheless, as pointed out in Q, when we use 
the effective coupling constant to describe hadronization corrections, it is not uniquely 
determined whether nic should be the quark or the hadron mass. As done for the purpose 
of the bottom-quark mass, we shall adopt a conservative choice and vary nic in the range 
1.5 GeV < rric < 2.1 GeV, that includes the current estimations for the charm pole mass 
as well as D-hadron masses |34|. Our default value will be mc=1.8 GeV. 



As for the DGLAP evolution operator, when evolving from fiop ~ rric to fip — Q, 
one typically crosses the bottom-quark mass threshold mi,, ref. [^] computed at NLO 
the matching conditions for the perturbative fragmentation function when crossing heavy- 
flavour thresholds. In our study, however, since we are working in the non-singlet ap- 
proximation and we are not accounting for gluon splitting and flavour mixing, we shall 
neglect such matching conditions. In fact, we checked that our results change very little 
according to whether we set in the DGLAP evolution operator, e.g, nj = 4ornj = 5as 
the number of active flavours. In any case, in our phenomenological analysis, whenever we 
have iiQF < mt, < hf, we shall implement the following factorized form for the non-singlet 
DGLAP evolution operator: 

En [as(MoF);«5'(/^F)] = ^^a^ ["sI/^of)' x En [asiml),as{i^F)] , (4.4) 

with = 4 and nj = 5 below and above the bottom-quark mass threshold, respectively. 
The 6-quark mass will be varied in the range 4.7 GeV < m^ < 5.3 GeV, as in with 
= 5 GeV being our default value. Elsewhere in our calculation, Uf will be consistently 
chosen according to the energy scale we are dealing with. 

4.1 Comparison with ALEPH data 

We shall first consider ALEPH data on D*^ production. As detailed in |31], such mesons 
can be in general produced from a ^ cc decay, from the decay of a primary 6-flavoured 
hadron produced in bb, from gluon splitting to cc or bb pairs, which subsequently 

hadronize or decay into a D*^ . The ALEPH Collaboration was able to subtract the 
Z^ — > bb and gluon-splitting contributions off and published the spectrum of D*^ mesons 
coming only from the cc primary source. In the following, we shall compare the predictions 
of our model with such a subsample, which will allow us to neglect secondary charm 
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production in the perturbative calculation as well as the singlet component of the DGLAP 
evolution operator, which in principle should play a role at LEP energies. Indeed, it was 



found out in |30| that implementing the singlet contribution does have an effect at xd < 0.4, 
and actually worsens the comparison with the ALEPH D*^ data coming from direct cc 
production (see figures 13 and 14 in ref. pO|]). 

In figure Q we present the spectrum given by our model, along with the D*^ ALEPH 
data, and investigate the dependence on the factorization scales fip and fiQp (figure ^ 
(a)), and on the choice of Q5(m^) and nic (figure |l] (b)). For the sake of comparison, 
both data and theoretical predictions are normalized to unity. As already observed in 
Q, the dependence on fioF, the scale entering in the initial condition of the perturbative 
fragmentation function, is fairly large, while the impact of the choice oi fip is pretty small. 
In particular, setting a lower value of fiQp, e.g. /^qf = "i,c/2, tends to deplete the small-xu 
region of the spectrum and to enhance the event fraction around the peak. Also, the peak 
is slightly shifted to higher xd if we choose fioF = mc/2. The prediction obtained for 
fiQP = 2mc reproduces quite well the Iow-xd data, while discrepancies are still present in 
the middle-high range. 

The dependence on Qs(m^) and rric is also quite relevant, as can be learned from 
figure m (b). In particular, a low value of rric, i.e. nic = 1.5 GeV, consistent with the quark 
mass rather than the Z)*^-meson mass, gives a pretty good description of the peak, but 
it worsens the comparison for xd > 0.7. On the contrary, a high value of rric, such as 2.1 
GeV, significantly moves the peak towards large xd and worsens the overall comparison. 
As for the effect of the variation of as{m'^), we find that it shifts the position of the peak: 
the lower as{m'^), the higher the value of xd at which the D spectrum is peaked. The 
dependence on the renormalization scales fiji and hqr is very little, and we do not present 
the corresponding plots for the sake of brevity. We also varied mb, the bottom-quark mass 



entering in eq. (4.4), but found out that it has negligible impact on the energy distribution. 



Overall, we can say that our model gives an acceptable description of the raise at low 
and average values of x^, while discrepancies are present around the peak, unless one sets 
a relatively low value for rric, and at very large xd- Our curves tend to be harder than 
the data and approach zero at large xd more rapidly. Although the comparison at very 
large xd is not completely satisfactory, using the In iZ-matching prescription, discussed 



in subsection 2.2., has nonetheless improved the spectrum near the endpoint xd = 1, 
as it is smoother and not oscillating any longer. We checked that if we had used the 
standard matching between NLO and resummed expressions as in Q, the charmed-meson 
distributions would have become negative for > 0.9. 

In any case, we are aware that our model, based on an extrapolation of perturbation 
theory, up to the replacement of the coupling constant as{k'^) asik"^), cannot be com- 
pletely reliable at very large xd- One can roughly estimate [0] x^j^max — 1 — ^/ruc — 0.85 
the maximum value of xd at which our model, or any model based on simple parametriza- 
tions of power corrections, such as the non-perturbative fragmentation functions ^, can 
be trusted. In fact, the authors of ref. [^0| managed to improve the comparison at large 
Xd, but they had to introduce a further free parameter which they tuned to data. In detail. 
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Figure 1: Comparison of the prediction yielded by our model with data from ALEPH on D*^ 
production a LEP. In (a) we investigate the dependence on pip (solid lines) and /ioF (dashes); in (b) 
the effect of the choice of as{m\) (solid) and rric (dashes). Such quantities are varied as discussed 
in the text. 



they replaced the Mellin variable N according to: 



1 + fN/N' ' ^ ' 

where A^' = exp[l/ (/3oas'(/U^))] in the coefficient function and N' = exp[l/ (2/3oa5(/^oi?))] 



in the initial condition, ref. |3C] used then / = 1.25 in its phenomenological analysis, 
as this choice led to good fits to charm-fragmentation data. In principle, we could also 
perform the replacement ( |4.5| ) and tune /. In fact, modifying the energy distribution 
at large xd will also have an impact at smaller values oi xq, since we have kept the 
normalization of data and theory curves fixed to unity. However, in this way our model 
would lose its crucial feature to be free from tunable parameters. Furthermore, given the 
theoretical uncertainty on our prediction, the value of / will depend on the particular set 
of perturbative parameters chosen. Therefore, for the time being, we prefer to stick to our 
parameter-free modelling of the hadronization and to postpone a more careful investigation 
of the very large-X£) regime of our spectra to future work, with the NNLO corrections 
implemented. In any case, we should never forget that, for the sake of consistency, whenever 
we modify the perturbative accuracy or the non-perturbative model, we should always 
reconsider the studies on B-hadron production and decay and check whether the results 
obtained in refs. Q and [Q] still hold. 

In the present analysis, as in ref. |^], we discard few points at very large xd and limit 
ourselves to xd < 0.85 when evaluating the from the comparison with the data. Even 
in this range, using our default values for the parameters in the parton-level computation, 
we are not able to acceptably reproduce the data, as we obtain x^/dof = 56.47/17. A 
better description of the data is nonetheless obtained if, e.g., we keep all quantities to their 
default values, but set hqf = "^mc (xVdof = 27.18/17) or as(m|) = 0.121 {x^/dof = 
30.52/17). Setting rUc = 1.5 GeV, we find xVdof = 32.29/17. As we are not fitting any 
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non-perturbative parameter to the data, such values of are acceptable. Also, they are 
of similar magnitude to those obtained in Q from the comparison with i?-hadron energy 
distributions at the pole for < 1 — A/nib- 

The overall impact of the inclusion of non-perturbative corrections at LEP energies via 
our model can be learned from figure ^, where we present our most significant predictions, 
i.e. the ones obtained with /ioF = 2mc and nic = 1-8 GeV (solid line), and with nop = 
nic and nic = 1.5 GeV (dotted), keeping the other quantities to their default values. 
In figure we also show the ALEPH D*~^ data and the purely perturbative results of 



ref. [11 1, where the authors used the standard coupling constant and resummed NLL soft 
and collinear contributions to both coefficient function and perturbative fragmentation 
function. The role played by power corrections is clearly remarkable throughout all x^- 
spectrum, and is essential to obtain an acceptable description of the data. In fact, the 



parton- level calculation of |11], which is the same as the one employed in |30|, needs to 
be convoluted with a non-perturbative fragmentation function to reproduce the data. We 
can also note in figure || that, while setting fioF = 'i.nic and rric = 1.8 GeV gives the lowest 
X^, the data around the peak are better described if we instead choose = 1.5 GeV and 
//OF = mc- 

Before closing this subsection, we remind that the possible reasons determining the 
fairly large theoretical uncertainties were already listed and detailed in 0. In particular, 
we have resummed large-x contributions to the coefficient function and initial condition in 
the NNLL approximation, but we have still matched the resummation to the NLO exact 
results, thus generating a mismatch between the NNLL terms in the resummed exponents 
(~ IxiN , etc.) and the remainder functions. We believe that the uncertainties should be 
milder if we used the exact NNLO results []l2| , p^ , 14, |l^. Moreover, lower theoretical errors 
should be expected if we also employed NNLL non-singlet DGLAP evolution equations, 
using NNLO non-singlet splitting functions |[T5[| . 

With respect to the analysis on i3-hadron production, the effect of the choice of scales 
and masses is here even more relevant: the dependence on such quantities is typically 
logarithmic, hence larger once they vary around rric rather than mi,. It is however interesting 
to notice that, unlike the comparison with the i?-hadron data, where setting //qf = i^T-b/'^ 
gave the best description of the data [Q], the charm- fragmentation data seem to prefer a 
quite high value of /UqFi since ^qf = 2mc yields the lowest x^- We believe that a full 
NNLO/NNLL analysis should clarify this issue as well. 

4.2 Comparison with CLEO and BELLE data 

We would like to compare the predictions of our model with the data on D^^ D*^ and 
D*^ production from the experiments CLEO [32| and BELLE collected at the T(45') 



resonance. In fact, since the value of the hard scale is much smaller than at LEP, such a 
comparison will help to shed light on the performance of our model and calculation at lower 



energies. Furthermore, ref. |30], using a NLO/NLL calculation and a non-perturbative 
fragmentation function with three parameters, managed to fit all i?-factory data, whereas 
some discrepancies were found with respect to the ALEPH data after evolving to LEP 
energies. Our case is clearly different, as our non-perturbative model is not tunable to 
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Figure 2: Results on charmed-hadron production at LEP {j — D), setting fiop = 2mc and 
771c — 1-8 GcV (solid), rric = 1.5 GeV and /ioF = Wc (dots), compared with the perturbative 
parton- level calculation of [0| (j — c, dashes) and the ALEPH D*^ data. 



data, but it will be nonetheless cumbersome to investigate how our predictions fare with 
respect to the different data sets at the T(4S') resonance and estimate the theoretical 
uncertainty. 

In figure ^ we present the comparison with CLEO and BELLE data on production, 
corrected for ISR effects. As pointed out in |Q, at the T{4S) resonance the contribution 
of cc pair production via gluon splitting is negligible, hence it is safe sticking to the non- 
singlet approximation of the DGLAP evolution equations, as done for analysis at LEP 
energies. The data sets which we consider are separately normalized to 1, for the sake of 
a consistent comparison with the theory curves, whose first moment reads, by definition, 
aN=i = 1. We vary renormalization and factorization scales, nic and asirn^) along the 
lines of our comparison with ALEPH. figure |3| (a) exhibits the dependence on fip and /ioFi 
figure ^ (b) the one on as{m?z) ^c- We do not present the effect of changing /xoi? 
and rrih, since it is very little, as already found at the pole. 

Unlike the comparison with the ALEPH data, where, though within the experimental 
and theoretical uncertainties, we succeeded in getting a reasonable fit of the data, our 
prediction lies quite far from the CLEO and BELLE data and there is no choice of 
parameters and scales, within our ranges, which can accommodate the experimental data. 
In fact, such data exhibit very small errors and, even if we limit our analysis to xd < 0.85, 
as we did before, we still obtain quite large typically x^/dof > 0(10). It is nonetheless 
interesting to notice that the best comparison is obtained for nic = 1-5 GeV: in this case, 
one is able at least to reproduce the rise of the spectrum up to x/j ~ 0.6, but still uncapable 
of describing the peak and the large- a; tail. As pointed out when comparing with ALEPH, 
a full NNLO/NNLL analysis is mandatory to reduce the theoretical error and should shed 
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light on the dependence on the quark (meson) mass as well. Ref. also presented 
data from CLEO and BELLE; the comparison with our predictions is however qualitatively 
similar to the one presented in figure |3| and we do not show it for brevity. 
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Figure 3: As in figure |l|, but comparing with _D" data from CLEO and BELLE experiments. In 
(a) we vary the factorization scales (solid) and ^qf (dashes); in (b) as(m^) (solid) and iric 
(dashes) . 

We present in figure |^ the comparison of the predictions yielded by our model with 
CLEO and BELLE data on Z)*''-meson production. The comparison with the CLEO D*^ 
data, which are affected by pretty small errors, is quite unsatisfactory and the values 
high. The BELLE spectrum instead exhibits larger errors, so that we are able to compare 
with the data at xd < 0.85 with quite small values. With our default parametrization, 
we obtain x^/dof = 45.23/36, while an even lower result, x^/dof = 32.10/36, is obtained 
if we set /Uqf = 2m-c, the same choice leading to the best fit to the ALEPH D*~^ spectrum. 

We finally show in figure I the comparison of our predictions with the data on D* 
production at CLEO and BELLE. Following [^], as far as the BELLE data are concerned, 
we present separately the spectra of the mesons decaying according to D*~^ and 
D*~^ — > D^, with the former presenting larger errors. The comparison is qualitatively 
similar to figure ^ with our model capable of describing well the data up to — 0.6, but 
failing to reproduce the peak and the large-xu tail. Drawing a parallel between figure ^ 
and figure |^, where our model, though within the uncertainties, led to a better comparison 
with respect to the ALEPH D*~^ data, one may argue that some major problems with our 
approach seem to appear once the process hard scale decreases. Later on, in section 6, 
we shall comment more about possible extensions of our calculation and non-perturbative 
model, which may eventually improve the comparison with the i?-factory data in x-space. 

As done when comparing with ALEPH, we show in figure ^ the prediction leading 
to the best fit to the S-factory data, i.e. the one obtained for fiQp = 2mc, along with 
the BELLE D*^ spectrum and the NLO/NLL perturbative prediction from ref. [^]. We 
note that the parton-level result is sharply peaked at large x, even more than in figure |2|: 
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Figure 4: Comparison of the prediction yielded by our effective-coupling model with D*^ data 
from BELLE and CLEO. The scales are varied as in figures and |. 




Figure 5: As in figures but comparing with the CLEO and BELLE data on D*^ production. 
When presenting the BELLE data, we distinguish the D*^ from the D*^ decay mode. 



in fact, the smaller phase space available at the T(45) resonance with respect to the 
pole enhances the probability of producing cc pairs near the threshold x = \. Overall, the 
impact of non-perturbative corrections in the coupling constant at the T(45') resonance 
looks even more important than at LEP energies. 

5. Results in A^-space 

In this section we present our analysis in Mellin space and compare our results with the 
experimental moments of the D-hadron cross section, measured by ALEPH, CLEO and 
BELLE. The data which we consider are the same as the ones which were analysed in the 



previous section in x-space. It was advocated in refs. |24, ^] that working in A^-space can 



be theoretically preferable, as one does not need any explicit form for the non-perturbative 
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Figure 6: The solid line is our best prediction for charmed-hadron production at the T(45') 
resonance (j — D); the dashed line is the purely perturbative c-quark spectrum (j — c) yielded by 
the computation in Also presented are the BELLE data on the D*° spectrum. 



fragmentation function and its moments can be fitted directly from the data. Of course, 
this issue does not apply to our case, since we are not fitting any parameter, but nonetheless 
it is still worthwhile to compare our results with the experimental moments. The A^-space 
investigation will be particularly interesting at the T(4S') resonance, where the x-space 
analysis has exhibited quite serious discrepancies. 

We calculate the moments of the Z)-hadron cross section directly from the A^-space 
formulas, i.e. eq. ( [4. 2D , and vary masses and scales as x-space, for the sake of estimating 
the theoretical uncertainty. The experimental moments are the same as the ones presented 
in ref. |3C], with the effect of electromagnetic initial-state radiation subtracted off. We just 
rescale them so that the first moment of all data sets reads o"7v=i = 1, as happens for our 
theoretical results. Following [^, we first evaluate the uncertainties on the moments due 
to to the variation of ^p, hqf, rric and as{m'^) separately, and then estimate the overall 
theoretical error summing in quadrature all individual uncertainties (see table 2 in ref. Q). 

The results of the comparison with the experimental moments are finally presented 
in figures ^[-|lO|, where we investigate how the prediction yielded by our model fares with 
respect to the moments of ALEPH £»*+ (figure 0), CLEO and BELLE (figure D, 
D*^ (figure |^ and D*~^ (figure ^) data. Since our model does not distinguish spin and 
electric charge, we shall always have the same theoretical moments, regardless of the kind 
of mesons we are comparing with. As found out in the x-space analysis, our predictions are 
affected by fairly large uncertainties; it is nonetheless interesting that, within the errors, 
the moments obtained using the resummed calculation provided with the effective-coupling 
model are compatible with the experimental ones. This result is especially remarkable for 
the comparison with the data collected at the T(4S') resonance (see figures |-|10|), which 
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Figure 7: Moments of charmed-hadron cross section according to our effective-coupling model 
(denoted by 'Theory'), compared with the moments of D*^ production at ALEPH. The theoretical 
errors are estimated by varying the parameters entering in the perturbative calculation, as discussed 
throughout the text. 
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Figure 8: As in figure ^ but comparing our prediction with the moments of the production 
cross section, measured by the BELLE and CLEO collaborations. 



exhibited instead relevant discrepancies in x-space. In fact, considering, e.g., the x-space 
spectra in figure |^ our model tends to underestimate the event fraction at small and 
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Figure 9: As in figures ^ and ^, witli our results faring against tlie iV-space D*° data from CLEO 
and BELLE. 
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Figure 10: Comparison of our tlieoretical prediction witli the iV-space D*^ data from CLEO and 
BELLE. As in the a;-space analysis, we plot separately the D*^ moments at BELLE, according to 
whether they decay via D*'^ — > D° or D*^ £>+. 



very large xo and overestimate the differential cross section for 0.6 < X£) < 0.8. Therefore, 
when evaluating integrated quantities like the moments, such effects get compensated and 
one is able to obtain a reasonable description of all A^-space data (figure ^). A similar result 
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was indeed found in |36], where parton showers and resummed calculations were used to 
describe i?-hadron production at the Z'^ pole. The spectrum yielded by the HERWIG 
Monte Carlo generator gave rise to a quite large x^/dof when comparing with h- 
fragmentation data in x-space, even after tuning a few parameters. However, HERWIG 
was able to reproduce fairly well the first few experimental moments of the B cross section. 

Referring, e.g., to the average value (x/j) of the spectra at 5-factories, correspond- 
ing to c7Ar=2, the experimental data yield (xd) = 0.610 ± 0.005 (CLEO) and 0.612 ± 0.006 
(BELLE). Using our default perturbative parametrization and effective-coupling model, we 
obtain instead (2;D)th = 0.647, above the experimental result. Different choices of the pa- 
rameters lead nevertheless to even larger or lower estimates of cfn=2- For example, setting 
rric = 2.1 GeV and the other parameters to their default values, we obtain (x£))th = 0.686, 
while for rric = 1.5 GeV the second moment reads (x£))th = 0.599, below the CLEO and 
BELLE data. In any case, as already pointed out for the purpose of the x-space analysis, 
a complete NNLO/NNLL should possibly decrease the theoretical error in Mellin space as 
well. Moreover, since even the best fits in x-space were obtained discarding the data points 
at X£) > 1 — A/ttt-c, we expect that the comparison with the experimental moments should 
eventually get worse for very large values of N , dominated by the x/j > 0.85 region in x- 
space, where our predictions are systematically below the experimental data. A hint for 
such a behaviour can be learned from figures |9| and IC, where already the (Jn=io theoretical 



moment is only marginally consistent with the data, even within the uncertainties. 
6. Conclusions 

We studied charm-quark fragmentation in e^e^ annihilation and used a recently pro- 
posed model, based on an effective strong coupling constant, as the only source of non- 
perturbative effects. Such a model was already employed in Q and gave a reasonable 
description of 6-quark fragmentation in e^e~ annihilation and some i?-meson decay data. 
We described charm-quark perturbative production following the perturbative fragmen- 
tation approach, with NLO coefficient function and initial condition of the perturbative 
fragmentation function, NLL DGLAP non-singlet evolution and NNLL large-x resumma- 
tion. Resummed expressions were matched to the exact NLO ones using the so-called 
In i2-prescription, which turned out to significantly improve the spectrum near the x = 1 
endpoint. The effective coupling was implemented in the NNLO approximation, as in 
refs. |,|]. 

We compared the predictions of our model with data from ALEPH, BELLE and CLEO, 
corrected for initial-state photon-radiation effects as in |30|. Throughout our analysis. 



since our non-perturbative model has no tunable parameter, we varied the quantities in 



the perturbative calculation within typical ranges, according to the values quoted in |34|| . 

We found that our model is able to acceptably describe, for x/j < 1 — A/rric and within 
the theoretical and experimental errors, the D*^ spectrum from ALEPH. In particular, 
the best fits to the data are obtained, within our chosen ranges, if we set the factorization 
scale entering in the initial condition to ^qf = 2mc. A value of rric consistent with the 
charm pole mass, rather than the D-meson mass, improves the comparison at small xd 
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and around the peak. Significant discrepancies were instead found with respect to the 
D^, D*^ and D*~^ data from i?-factories, where we succeeded in obtaining x^/dof ~ 1 
only when comparing with the BELLE D*^ spectrum, affected by pretty large errors. The 
experimental data on and D*^ production at CLEO and BELLE and on D*^ at CLEO 
exhibit instead very small errors and we did not manage to obtain a reasonable x^/dof, 
even within the theory error. We just noticed that setting = 1.5 GeV gives a good 
description of the data for xd < 0.6, but some major disagreement is still present for larger 
values oi Xd- In Mellin space, however, within the fairly large theoretical uncertainties, we 
managed to reproduce the first ten moments of all considered data samples. We expect 
nonetheless that for larger values of N some discrepancy should appear, consistently with 
the observation that even the best rc-space fits were obtained discarding few large- data 
points. 

Anyway, we remind that some problems with reproducing both ALEPH and i3-factory 
data were already encountered in ref. |30|, where the authors employed a NLO/NLL calcu- 
lation, a non-perturbative fragmentation function with three free parameters and rescaled 
N according to eq. ([4.5|). The hadronization model was tuned to reproduce fairly well all 
S-factory data on neutral as well as charged D- and L)*-meson production. Nevertheless, 
the best-fit parametrization did not succeed in reproducing the ALEPH D*^ spectrum 
after evolving to LEP energies. The conclusion of the analysis carried out in [^] was 
that, in order to reconcile both LEP and i?-factory data, it was necessary to include power 
corrections in the process-dependent coefficient function, depending on the process hard 
scale. This way, one should be able to describe all data, still using the same perturbative 
accuracy and the same functional form for the non-perturbative part. However, due to the 
errors in the intermediated xd region, ref. [^0| was not able to discriminate whether the 
missing power corrections should behave according to a 1/Q or a l/Q^ power law, with 
Q being the centre-of mass energy. In any case, being mz much larger than ?tt.x(45'), one 
should expect that such a power correction should mainly modify the spectra at the T(45) 
resonance, in such a way that the fits of the non-perturbative fragmentation functions at 
the pole, presented in |30|, should eventually work even at 5-factory energies, after 
minimal adjustments |]38| . 

As far as our work is concerned, we do find it interesting that, although within the 
theoretical and experimental errors and after discarding few data points at very large x^, 
our parameter- free model yields x^/dof ~ 1 from the comparison with ALEPH D*^ data 
and reproduces the moments of all analysed data sets. The discrepancies of our prediction 
with respect to the very precise data from CLEO and BELLE in x-space clearly deserve 
further investigation. The results in this paper, along with the ones reported in Q, seem 
to indicate that the model works better for heavy-quark fragmentation at the pole, 
while more serious problems show up once the hard scale is lowered. However, only a 
power correction, such as the one understood in [^], mostly relevant at large x or A^, may 
not to be enough to solve the discrepancy with the i3-factory data, as figures ^-js] show 
disagreement even around the peak and at small xd- 

The theoretical uncertainty is expected to decrease after the inclusion of NNLO coef- 
ficient functions initial condition |14, 16| and non-singlet splitting functions |15|, 
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which win also promote DGLAP evolution to NNLL accuracy in the non-singlet sector. 
Since within our approach we are including power corrections in an effective coupling, any 
perturbative improvement, such as accounting for 0(a|) contributions, will necessarily im- 
ply the inclusion of non-perturbative corrections as well. Furthermore, the change ( 3.12| ) in 
the coefficient A^'^^ has been implemented in the threshold NNLL expressions, but not yet 
in the splitting functions, whose NNLO corrections do contain a contribution ~ ^4^^^ ]l5| . 
Including such a term in the splitting functions, along with the redefinition A^^^ — > A^^\ 
may shift the xd spectrum and possibly improve the comparison with the i?-factory data, 
as found in Q for 5-hadron production at LEP and SLD. We should also expect a relevant 
impact on our analysis of the possible inclusion of large-x next-to-next-to-next-to-leading 
logarithmic (NNNLL) terms, whose coefficients have been denoted by A^'^\ B^^) ^(3) 
in Eqs. ( |2.7D and ( [2.111) . All such coefficients will be modified in a fashion analogous to 
eq. when using asik"^)- Moreover, the implementation of higher-order threshold 

contributions in the exponents ( |2.7| ) and (2.11) will also lead to the inclusion of further 
power corrections since, as part of our model, we performed the Mellin transforms exactly. 

The other guideline to obtain better agreement with the x-space data consists in 
modifying the effective-coupling model, e.g., introducing a correcting term 6ds{k'^) as in 
eq. ( tj.l4 ), possibly containing extra parameters. In order to accommodate both D and 
D* data, 6as{k^) may possibly depend on the spin of the considered hadron. However, 
before speculating about its functional form, we believe that we still need a NNLO/NNLL 
calculation to reduce the theoretical uncertainty and deal with a more stable prediction. In 
fact, without a NNLO/NNLL analysis, function dasik"^) will largely depend on the values 
chosen for the perturbative parameters and considerably vary according to whether, e.g., 
one sets rric = 1.5 or 2.1 GeV, //qf = "ic/2 or 2mc, and so on. 

The laige-XD behaviour of our spectra may be improved as well, since even the best 
comparisons with the data were obtained in this paper for < 1 — A/tti-c and in Q for 

< 1 — A/mfe. An option could be the prescription ( [4.5| ) suggested in |Q; in fact, any 
modification at large xd will indirectly affect, via normalization, the energy distribution 
at smaller as well. Nevertheless, once again, given the uncertainties exhibited by 
our predictions even at large xd, this investigation should be better performed using a 
calculation of higher accuracy. 

Ideally, once the above issues are clarified, one may think of using our model to describe 
D- and 5-hadron production at the Tevatron accelerator, along the lines of refs. p9| , 40|, 
and extend the results to LHC energies. Nevertheless, unlike the standard analyses, where 
a non-perturbative fragmentation function is fitted to e~^e~ data and then used in the 
hadron-collider environment, we are not tuning any parameter to the e~^e~ data. Therefore, 
possible studies at hadron colliders will be independent checks of the capability of our model 
to reproduce heavy-quark fragmentation data. 

Moreover, we can use the NLO perturbative calculations in | |4l| , along with the 
effective coupling constant, to predict bottomed-hadron spectra in top (i hW) or Higgs 
[H — > hb) decays at the Tevatron and LHC. Finally, the c- fragmentation results here pre- 
sented can be compared with the ones yielded by Monte Carlo generators, extending the 
analysis carried out in [EgI, where parton shower algorithms and resummations were com- 
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pared for the purpose of i?-hadron production in e^e annihilation, top and Higgs decays. 



For such a comparison to be consistent, however, even the HERWIG [^] and PYTHIA |43| 
generators will have to be tuned to the same LEP and S-factory data analysed throughout 
this paper. It will also be very interesting to implement the effective coupling constant 
to replace, e.g., the cluster model |44] which simulates the hadronization in HERWIG and 
investigate how the Monte Carlo results fare with respect to the experimental data on D- 
and i?-hadron production. This is in progress as well. 
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